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Abstract 

In this paper we examine the possibility of testing the equivalence principle, in its weak form, by analyzing the 
orbital motion of a pair of artificial satellites of different composition moving along orbits of identical shape 
and size in the gravitational field of Earth. It turns out that the obtainable level of accuracy is, realistically, 
of the order of 10^^° or slightly better. It is limited mainly by the fact that, due to the imavoidable orbital 
injection errors, it would not be possible to insert the satellites in orbits with exactly the same radius and that 
such difference could be known only with a finite precision. The present-day level of accuracy, obtained with 
torsion balance Earth-based measurements and the analysis of Earth-Moon motion in the gravitational field 
of Sun with the Lunar Laser Ranging technique, is of the order of 10^^^. The proposed space-based missions 
STEP, ^SCOPE, GG and SEE aim to reach a lO^^^-lO-i** precision level. 
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1 Introduction 

The Weak Equivalence Principle (WEP) and the Einstein Equivalence Prnciple (EEP) (Will 
1993; Ciufolini and Wheeler 1995; Haugan and Lammerzahl 2001; Will 2001) are the corner- 
stones of Einstein General Theory of Relativity (GTR) and of all the other competing metric 
theories of gravity. The WEP states that non-rotating, uncharged bodies of different structure 
and compositions and with negligible amount of gravitational binding energy per unit mass fall 
with the same acceleration in a given gravitational field, provided that no other forces act on 
them. The EEP (called also Medium Strong Equivalence Principle) states that the outcome of 
any local, non-gravitational test experiment is independent of where and when in the gravita- 
tional field the experiment is performed. If the gravitational binding energy per unit mass of 
the freely falling bodies in a given external gravitational field is not negligible, as is the case for 
astronomical bodies, we have the Strong Equivalence Priciple (SEP. It is called also Very Strong 
Equivalence Principle). When the self-gravity of the falling bodies is accounted for it might 
happen, in principle, that their accelerations are different. This violation of the universality of 
free fall, called Nordvedt effect (Nordvedt 1968a), is predicted by all metric theories of gravity 
apart from GTR which, instead, satisfies the SEP. Lunar Laser Ranging (LLR) is able to test 
this effect for the motion of Earth and Moon in the gravitational field of Sun (Nordvedt 1968b). 
No different accelerations have been found at a 10~^^ level (Anderson and Williams 2001). Of 
course, this is also a test of the WEP. 

To go to space, where much larger distances, velocities, gravitational potential differences 
with respect to Earth, and, most important, free fall for an, in principle, infinitely long time are 
available, is important to perform very accurate tests of post-Newtonian gravity. In this paper 
we wish to investigate the level of accuracy which could be reached in testing the WEP, by 
analysing the orbital motion of a pair of Earth artificial satellites of different compositions for 
which the effects of self-gravity would be, of course, negligible: it is analogous to the analysis of 
the motion of Earth and Moon in the gravitational field of Sun in order to test the SEP. For the 
already performed experimental tests of the WEP on Earth see (Will 2001) and the references 
therein. The present-day level of accuracy is 4 x 10~^^ (Bassler et al 1999). It has been obtained 
in the experiments of the so-called Eot-Wash group of the Washington University by means of 
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a sophisticated torsion balance in an Earth based laboratory set-up. Other proposed space- 
based experiments are the very complex and expensive STEP (Lockerbie et al 2001), /^SCOPE 
(Touboul 2001), GG (Nobih et al 2000) and SEE (Sanders et al 2000) missions^ whose goal is 
to reach the 10~^^-10"^^ level of accuracy. 

The reason for searching for more and more accurate tests of the equivalence principle resides 
in the fact that all approaches to quantizing gravity and to unifying it with the other funda- 
mental interactions currently under study are capable of predicting violations of the equivalence 
principle at some level. For example, departures from universal free fall accelerations of the 
order of 10~^^ have been calculated in (Damour and Polyakov 1994a; 1994b) in the string theory 
context. Violations of the WEP are predicted also by nonsymmetric theories of gravity (Will 
1989). In (Moffat and Gillies 2002) a violation of the WEP of the order of 10"^^ is predicted 
in the context of a non-local quantum gravity theory. 

2 The orbital period 

In this section we examine the possibility of testing the WEP by measuring the orbital periods 
of a pair of Earth artificial satellites of different compositions. It can be thought of as a 
comparison between two pendulums with enormously long threads swinging for an extremely 
long time^. 

The orbital period of a satellite freely falling in Earth's gravitational field can be written as 

2vr 27r, , , 

T=— — l-An, (1) 

n(l + An) ^ ' 

in which the Keplerian unperturbed period is 



""^See also on the WEB http://einstein.stanford.edu/STEP/, http://www.oncra.fr/niicroscope/ 
Ihttp:/ /tycho. dm. unipi.it/~nobih/ggproject. html and http://www.phys.utk.edu/scc/ Notice that while 
the first three missions are in advanced stages of planning and hardware testing, and are expected to be 
launched in the next few years, SEE is still undergoing rigorous conceptual evaluation and is not yet a scheduled 
mission. 

^In the case of the SEP-LLR experiment the amplitude of the parallactic inequality long-period harmonic 
perturbation, proportional to cosZ?, where D is the synodic phase from New Moon, is sensitive to the possible 
different falling rates of Earth and Moon toward the Sun. 
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where G, M and a are the Newtonian gravitational constant, the mass of Earth and the 
satelhte semimajor axis, respectively. We have explicitly written the square root of the ratio of 
the incrtial to the (passive) gravitational mass of satellite The quantity An represents the 
various kind of perturbations, of gravitational and non-gravitational origin, which affect n. For 
example, the even zonal harmonic coefficients of the multipolar expansion of Earth gravitational 
potential, called geopotential, induce secular perturbations on n. The most important one is 
that due to the first even zonal harmonic J2 and, for a circular orbit with eccentricity e = 0, it 
is given by 

Anilf = -7^2(-V(l-3cos^z), (3) 



a / 

where R is the Earth radius and i is the inclination of the orbital plane to the equator. Also 
the time-varying part of Earth gravitational field should be considered, in principle, because 
the Earth solid and ocean tides (lorio 2001a) induce long-period harmonic perturbations on n. 
For a given tidal line of frequency / the perturbations induced by the solid tides, which are the 
most effective in affecting the orbits of a satellite, can be written as 



m=0 \ -n- / \ " / p=o q=-oo 



e^) dGipq 



2{e+l)Gepq- cos -ifimpq , 

(4) 

where Hp are the tidal heights, kf^ are the Love numbers. Aim — ^ ''^^^(i+tn 'yr^' ^irnp{i) and 
Gipq{e) are the inclination and the eccentricity functions (Kaula 1966), respectively, and '-fjempq 
is the frequency of the tidal perturbation built up with the frequencies of the lunisolar variables 
and of the satellite's orbital elements. 

In order to test the WEP, we propose to measure, after many revolutions, the difference 
of time spans which are multiple N of the orbital periods ATjv = Tj^^ — Tj^^ of a couple of 
satellites of different composition orbiting the Earth along circular orbits of almost same radius 
a 

a^^^ = a, a^^) = a + d. (5) 

The small difference d, which, in principle, should be equal to zero, is due to the unavoidable 
orbital injection errors. It can be made extremely small with a rocket launcher of good quality: 
for example, at the beginning of their mission the semimajor axes of the two GRACE spacecrafts 
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were different by an amount of just'^ 0.5 km. 

The observable quantity in which we are interested is Aip to be measured after orbital 
revolutions, where 

i)^^^ = 11) = !, ^^^^ = ^ + ^11), (6) 

with Alp /if) -C 1. The quantity Aip can be expressed in terms of the standard Eotvos parameter 
7]. Indeed, the inertial mass of a body is composed by many types of mass-energy: rest energy, 
electromagnetic energy, weak-interaction energy, and so on. If one of these forms of energy 
contributes to differently than it does to m; we can put (Will 1993) 



m„ 



m 



(7) 



where is the internal energy of the body generated by the interaction A and t]^ is a di- 
mensionless parameter that measures the strength of the violation of WEP induced by that 
interaction. Then 
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From eq.® it can be obtained 
so that, from eq.® 
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From eqs.dH)-® we can write 

ATn = N2n 



An 



(1) 

obi 



An 



{2)' 
obi 



n 



(2) 



n 



(1) 



n 



(1) 



n 



(2) 



(9) 
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In eq. fllH) we consider only the gravitational even zonal perturbations due to Earth oblateness, 
as usually done in orbital reduction programs in which a reference orbit including the J2 effects 
is adopted^. They can be summarized as 



An, 



obi 



a 



(12) 



•^See on the WEB "http://www.csr.utexas.edu/grace/newsletter/2002/augtist2002.htm] 
^This approximation will be justified later. 
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where the Qi^ = Qe{h e; Je) functions include the even zonal harmonics coefficients Ji, the 
eccentricity e, the inclination angle i and some numerical constants. For example, for i = 2 
and e = 0, Q2 = — 1</2(1 — 3 cos^ i). By using the expansions of eqs.®-® it is possible to solve 
eq.jni) for^ AtpN 

^-ATn -{A + B) 



Alp, 



N 



/GM 
Ntt 



{A + C) + {B + D) 



with 



( l-2(?^ 



(13) 

(14) 
(15) 

(16) 
(17) 

Eq.lfT^. together with eqs.(fT ^ -(fT7 |) . allows to evaluate the accuracy obtainable in measuring 
the quantity Aip]^. The error in measuring the difference of multiples of the orbital periods 
yields 



A 


= 3d^/a, 




B 




r, / 1 + 21 

3£d^a-y — 




1=2 




C 


= 2v^, 




D 







[5iAM] 



VGM 6{ATj 



Nn {A + C) + {B + D)' ' ' 

while the error in Earth GM, which amounts to 8x10^^ cm^ s~^ (McCarthy 1996), yields 

= N2.VGM{A + C) + {B + Dy 

The uncertainty in Earth even zonal harmonics Ji, of which 6J2 amounts to 7.9626x10^^^ 
(Lemoine et al 1998), has an impact given by 

-If X p + C) + (5 + D)] + (A + 5) X ^ 



[{A + C) + {B + D)]^ 
The errors in AipN due to the uncertainties in a and in d are 

X [{A + C) + {B + D)] + {A + B)x d\iA+c)HB+D)] 



X 6Je 



(20) 



[S{A^n)]c 



^Here we have considered i'-^-' 



,(2) 



[iA + C) + iB + D)]^ 

Moreover, notice that if the two sateUites would be inserted in 



X 5a,(21) 



counter-rotating orbits, in ATjv it must be included also the time shift AT^™ oc TV x 47r^^^ due to the general 
relativistic gravitomagnetic clock effect (Mashhoon et al 1999; 2001; lorio et al 2002) induced by the off-diagonal 
components of the metric proportional to the proper angular momentum J of Earth. Indeed, it turns out (lorio 
et al 2002) that such time shift is independent of the ratio of the inertial to the gravitational masses of the 
satellite. It would act as a lower limit of the order of x lO"'' s to (Ar;v)cxp if it was detectable. 
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[S{AM]d 



-^x(C + D) 
[{A + C) + {B + D)f 



X Sd. 



(22) 



Notice that [(5(A'0jv)] j^, [S{AipN)]a and [(5(A'0jv)]d do not depend explicitly on the number 
of orbital revolutions N. In order to calculate them we need the explicit expressions of the 
derivatives with respect to 
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and those with respect to d 
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(29) 
(30) 



In order to fix the ideas, let us consider the orbit of the proposed LARES laser-ranged 
satellite with a — 1221 Q km, i—lQ deg. Let us assume (Peterson 1997) o? = 5 km. With these 
data we have, for ^ = 2, 

<5(AT^) 



N ATj 



N 



7.3 X 10"^ s"^ X 



-14 „-l 



[(5(A^^)]gm = 7.4 X 10-- s ^ X 



N ' 
(AT^) 

exp 



AT 



[5(AV'iv)]a = 5 X 10-^^ cm-^ X 5a, 
[5{AiljN)]d = 1 X 10"*^ cm-^ X 5d. 



(31) 

(32) 
(33) 
(34) 

(35) 
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The accuracy in measuring the difference of the muhiples of the orbital periods of the satelhtes is 
a crucial factor in obtaining a high precision in Aip. However, it could be possible, in principle, 
to choose an observational time span covering a very high number of orbital revolutions. It 
should be noted that 5(AT/v) accounts for both the measurement errors and the systematical 
errors induced by various gravitational and non-gravitational aliasing phenomena. The latter 
ones play a very important role in strongly limiting the possibility of measuring, e.g., the 
gravitomagnetic clock effect: it amounts to 10""^ s after one orbital revolution for orbits with 
e = z = while in (lorio 2001b) it turned out that the systematic errors induced by the present- 
day level of knowledge of the terrestrial gravitational field are up to 2-3 order of magnitude 
larger. So, it should not be unrealistic to consider (5(A77v)sys ~ 

10"^ - 10"^ s. This would 

imply [5(A^^)]at, ~ icr!^_ 

The error due to the difference d in the semimajor axes of the satellites, which turns out to 
be the major limiting factor, cannot be reduced, in principle, by waiting for a sufficiently high 
number of orbital revolutions because it is independent of A^. A small improvement could be 
obtained with the use of a larger semimajor axis. A geostationary orbit with a = 42160 km 
would allow to get [5(AV'7v)]at^. = 1 x 10"^ s'^ x and [(5(AV^;v)]d = 3 x lO'^o cm'^ x 5d. 

However, in this case, for a fixed time span, we would have at our disposal a smaller A^. In 
regard to the systematic part of the error 6d, it should be noted that there are no secular 
or long-period perturbations of gravitational origin on the semimajor axis of a satellite. The 
non-gravitational perturbations could be reduced to a good level by adopting the drag-free 
technology. In regard to Sd^^p, it is important to note that in the present GRACE mission 
(Davis et al 1999), making use of a K/Ka-band intersatellite link that provides dual one-way 
range measurements, changes in the distance of the two spacecrafts can be established with an 
accuracy of about 10~^ cm or even better. 

It is interesting to note that if we neglect in all calculations B and D and the related 
derivatives, i.e., if we neglect the effects of Earth oblateness, it turns out that the numerical 
results of eqs. ()3H) - (jn3j) do not change. It is very important because it means that our choice of 
neglecting the contribution of the non-gravitational perturbations in An is a posteriori correct. 
Indeed, the perturbing acceleration due to Earth J2 on LAGEOS is of the order of 10~^ cm s~^, 
while the impact of the direct solar radiation pressure, which is the largest non-gravitational 
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perturbation on LAGEOS, amounts to 10~^ cm s"^ (Milani et al 1987). Moreover, these 
conclusions imply that the errors in the inclination i, which enters the even zonal harmonic 
perturbations due to the geopotential and the non-gravitational perturbations, can be safely 
neglected, as done here. The same considerations hold also for the tidal perturbations: suffices 
it to say that the effects of the 18.6-year and the Ki tides, which are the most powerful 
in perturbing the satellite orbits, on a LAGEOS-type sateUite are six orders of magnitude 
smaller than those due to the static J2 even zonal part of geopotential. Also the tiny general 
relativistic gravitoelectric correction to the orbital period induced by the Schwarzschild part of 



the metric^, which depends on ymi/mg, can be neglected because for a LAGEOS-type satellite 
the disturbing acceleration is of the order of 9xl0~^ cm s~^ (Milani et al 1987). 

3 The longitude of the ascending node 

The longitude of the ascending node O is one of the best accurately measured Keplerian orbital 
elements of Earth artificial satellites. Then, we wish to examine if it would be possible to use 
it in order to test the equivalence principle. 

Let us recall that there are two kinds of long-period perturbations on the node of an Earth 
satellite. First, the static oblateness of Earth induces a secular precession of Vt through the 
even zonal harmonics of the geopotential. Second, the time-varying part of Earth gravitational 
potential induces tidal harmonic perturbations on O (lorio 2001a). Then, we can pose, by 
including the solid Earth tidal perturbations 



a 



R J p=oq=-oo sm^Vl-e^ 



(36) 



where, ior £ — 2 and e — 0, Q2 — — §J2C0si. By considering a couple of satellites of different 
compositions freely orbiting along almost identical orbits we could measure the difference of 
their secular nodal rates AO = f2^^^ — In this case, by posing ^ = ^J^, the violating 

parameter = expressed as 



An 



+ A + B 



AC = , (37) 



^For a circular orbit it is given by T(°) x Oge = ^ ^ttVCMo (Mashhoon et al 2001). 
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with^ 

A = ^i?^A6;,a-('^), (38) 



B = dY,R%( 

1=2 \ 



^)a-(-), (39) 

c = j:r'g^^'^''^i (40) 

(41) 



e=2 



where AQ^ = Qf^ — Qf^ is the difference in the Qi functions of the two satelhtes induced by 
the inclinations and the eccentricities. The derivatives with respect to J^, a and d are 

dA ^^.^aAft\ _(^)_ ^^^^ 



=2 




dB 

die 
dC_ 

die 

OA 

da 

dB f3 + 2i\ f5 + 2i\ _(7+2i 

da 
dC 
da 



^E^M^)(^)«-^^^ (43) 

(44) 

(45) 

(46) 

(47) 



For i = 2 and a GPS orbit, by assuming a = 26578 km, i = 55 deg and d = 5 km the errors in 
are^ 

mO]An = ^yj^c- B) ^^"" (mas/yr)-^ x 6{An), (49) 



^Here we neglect the non-gravitational perturbations on the nodes. Contrary to, e.g., the perigees us, the 
nodes are rather insensitive to such non-geodesic accelerations (Lucchesi 2001; 2002). In regard to the Earth 
solid tides, they have been neglected because their impact is several orders of magnitude smaller. For example, 
the amplitude of the nodal rate perturbation induced by the Ki tide is five orders of magnitude smaller than 
that due to the even zonal harmonic coefficient J2 of the geopotential for a GPS orbit. 

®In the calculations it turns out that the effect of AQ2 in A and for Ai — 1 deg, can be ne- 
glected. On the other hand, with a good quality rocket launcher it is possible to insert two spacecrafts 
in the same orbital planes up to lO"** deg, as in the case of the GRACE mission. Sec on the WEB 
|http://www. csr.utexas.edu/grace/newsletter/2002/august2002. html 
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exp 



2^{GMf{C - B) 



X 5{GM) - 3 X 10"'" (deg/day)"' x (AO) 



exp) 



^x(C-i3) 



d{C-B) 
d.J2 



{C - B f 

^x(C-ig)-^x(^ + g) 



X 5J2 = 2 X 10 
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(50) 
(51) 



f x(C + ^) 



(C - BY 
X 5d = 1 X 10~^ cm~^ X 5d. 



X 5a = 2 X 10-^^ cm-^ x 5aj(52) 

(53) 



(C - B)2 

It can be noticed that the major hmiting factor is the term due to the error in the difference 
of the nodal rates 5{AVt) = + SQ^^\ Indeed, the experimental error in measuring the 

secular rate of the node is of the order of 1 mas yr^^. The systematic error in Q due to the 
uncertainty on J2 is, for a GPS satellite, almost 3 mas yr~^. In the case of the error in d, the 
same considerations as for the orbital periods hold. 

4 Conclusions 

In this paper we have shown that the comparison of the orbital motions of a pair of artificial 
satellites of different compositions moving along identical orbits in the gravitational field of 
Earth in order to test the Weak Equivalence Principle is not competitive with the already 
performed tests with torsion balances on Earth and the Lunar Laser Ranging technique, and 
the dedicated space-based missions STEP, GG, /iSCOPE and SEE. 

We have considered the orbital periods and the secular nodal rates. The analysis of the 
orbital periods seems to yield more precise measurements. The major limiting factor is rep- 
resented by the difference in the orbital radiuses induced by the unavoidable orbital injection 
errors and the related uncertainty. By assuming 5d < 1 cm or less the achievable precision is 
of the order of 10-^°-10-^^ 
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